Screening in the two-dimensional electron gas with spin-orbit coupling 
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We study the interplay between electron-electron interaction and Rashba spin-orbit coupling in 
the two-dimensional electron gas. Using the random phase approximation we predict new screening 
properties of the electron gas which result in extension of the region of particle-hole excitations and 
in spin-orbit induced suppression of collective modes. In particular, we observe that the spin-orbit 
coupling gives a finite lifetime to plasmons and longitudinal optical phonons, which can be resolved 
in inelastic Raman scattering. We evaluate the experimentally measurable dynamic structure factor 
and estimate the range of parameters where the described phenomena are mostly pronounced. 
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Fundamental issues of interaction effects in the two- 
dimensional electron gas (2DEG) are at the center of dis- 
cussions since the early days of its fabrication . Quite 
generally, screening described by the dielectric function 
forms the basis for understanding a variety of static and 
dynamic many-body effects in electron systems Q. The 
dielectric function of the 2DEG was computed a long 
time ago by F. Stern [j| within the random phase ap- 
proximation (RPA) . Different quasiparticle and collective 
(plasma) properties deduced from these expressions were 
confirmed experimentally soon after Q. Recent experi- 
ments measuring plasmon dispersion, retardation effects 
and damping |a, la III unambiguously show the impor- 
tance of correlations between electrons. 

More recently, the possibility of manipulating spin in 
2DEG by nonmagnetic means has generated a lot of ac- 
tivity #j. The key ingredient is spin-orbit (SO) coupling 
tunable by an applied electric field [2( . A recent example 
for SO-induced phenomena is the spin-Hall effect jlCTI ] - 

In this context understanding of the combined ef- 
fects of electron-electron correlations and SO coupling in 
2DEG becomes vitally important. In the present paper 
we study the RPA dielectric function at finite momenta 
and frequencies treating the SO coupling exactly. We 
observe two interesting phenomena: an extension of the 
continuum of particle-hole excitations and SO-induced 
damping of collective modes. New features in the dy- 
namic structure factor are predicted, and SO-induced 
widths of both plasmon and longitudinal optical phonon 
spectra are evaluated. Their values lie in the ranges 
which are experimentally accessible by now in the inelas- 
tic Raman scattering measurements. We also revisit the 
derivation of zero-frequency limit of the dielectric func- 
tion and present an analytic result which does not contain 
an anomaly at small momenta 

We consider a 2DEG with SO coupling of the Rashba 
type \l2\ described by the single-particle Hamiltonian 
H = + ctR^cr x k), where n is a unit vector normal 
to the plane of 2DEG and h = 1. The dispersion relation 
is SO-split into two subbands e k = 5— ^ ± cvrA; with two 



distinct Fermi momenta k± = kp =p m*an = kp Ur, 
but with the same Fermi velocity vf = kp/m*. The 
effective Coulomb interaction is V^' = V q /s qu , where 
V q = 2-Ke 2 1 '(qeoo), £00 is the (high-frequency) dielectric 
constant of medium, and the dielectric function e quJ = 
describes effects of dynamic screening. In 
RPA, e qu} = 1 — VqTl quJ , where is the polarization op- 
erator. In the presence of SO coupling U qLJ — II+, + 11^ 
includes contributions from intersubband (+) and intra- 
subband (— ) transitions: 
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where the form factors f^ k+q = |[1 ± cos(</>k — </>k+ q )] 
originate from the rotation to the eigenvector basis. 

We now discuss our results for e^f A . We have derived 
them in the analytic form which will be presented in a 
subsequent paper. Here we discuss the main SO-induced 
effects for the range of small z and w, and describe how 
they depend on the SO-coupling strength y — kji/kp 
and on the Wigner-Seitz parameter r s = 1/ \Jrma*g 
which is defined by an electron density n and effec- 
tive Bohr radius a* B = EooTOeas/m*. Our main results 
are compactly summarized in Fig. ^ which shows the 
contour plots and cross sections of the structure factor 
S(z, w) — — Im [l/e(z, w)] for different values of y and r s . 
The introduced notations z = q/2kp and w — m*u>/2k F , 
are, respectively, the dimensionless wave vector and fre- 
quency. 

First of all, we observe an exstension of the contin- 
uum of particle-hole excitations (or Landau damping 
region) defined by £2 7^ 0. A new wedge-shaped re- 
gion of damping is generated by the SO coupling. It 
is bounded by two parabolas — (z — y) 2 — (z — y) < w < 
(z + y) 2 + (z + y) (black-dashed lines in Fig. ^) and at- 
tached to the boundary w = z 2 + z (white-dashed line in 
Fig. ^) of the conventional particle-hole continuum (ob- 
tained in the absence of SO coupling according to 
Another boundary w — z 2 — z of the latter transforms 
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FIG. 1: Contour plots of S(z,w) showing the SO-induced 
wedge-shaped damping region (bounded by the dashed lines). 
The plasmon mode is depicted by the bold line. Insets show 
the cross-sections S(z = 0.1, w) as a function of w. Left panel: 
y = 0.07, r s = 0.2. Right panel: y = 0.0A, r s = 0.6. 



into w — [z — y) 2 — (z — y) for nonzero y (this occurs at 
z > 1 which are not shown in Fig.QJ. The new regions of 
damping reflect an opened possibility for transitions be- 
tween SO-split subbands, and S(z, w) strongly depends 
on y inside these regions. Their boundaries can be de- 
termined from a simple consideration of extremes of the 
denominators in On the other hand, within the con- 
ventional boundaries w = z 2 ± z the function S(z, w) is 
modified by SO coupling only slightly and can be approx- 
imated by Stern's expressions [3(- 

The second SO-induced modification is the broadening 
of a plasmon collective mode within the new particle- 
hole continuum. We determine numerically the position 
of the plasmon spectrum from the equation ei(z,w) = 
and depict it in Fig. ^ by the bold line, in black where 
£2 = (undamped plasmon, the structure factor is delta- 
peaked) and in gray where £2 7^ (SO-damped plasmon, 
the structure factor has finite height and width). Two 
different cases are possible: (I) the plasmon enters only 
once into the SO-induced damping region (Fig. 2] left 
panel); (II) it enters twice escaping for a while after the 
first entrance (Fig. ^ right panel) . We observe that the 
position of the plasmon in the presence of SO coupling is 
almost indistinguishable (at least, on the scale of Fig. 
from the RPA result derived in the absence of SO cou- 
pling EH 



zjz + 2fg) Af 2 + Af s z 3 + 2 4 nf , 
w p i(z) = — x j — ; — N Q(z -z), (2) 



2f, 



z(z + 4f s ) 



where f s = r s A/8 and z* is the real positive root of 
the equation z 2 (z + Af s ) = Af 2 . Later we will comment 
more on how small might be the difference between the 
exact SO-modified plasmon dispersion and that given by 
(0). At the moment we will use in order to calculate 
critical values of y and r s which separate the cases (I) 
and (II). The result is represented in Fig. [3 where the 
plane (y, r s ) is divided into the respective domains I and 
II. Changing y and r s , we can tune the relative position 



of the new particle- hole region and w p i(z). 

In the new region of damping the dynamic structure 
factor is very well approximated by the Lorentzian 



S{Z, W) SO-damp pi 



a(z) 



j(z) 



(w-w p i(z)) 2 +7 2 (z)' 



(3) 



which describes the SO-damped plasmon with the width 
j(z). The weight factor a(z) — Tr(dex/dw\ w — w ^z))" 1 
is practically independent of y and given by its value at 
y = Q 



a{z) 



/I[16ff-z 4 (z + 4f s ) 2 



2f sy /(4f 2 + Af s z 3 + z A )(z + Af s ) 3 
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At the same time, the SO-induced width 'y(z) = 
a(z)e2(z,Wpi(z))/ir strongly depends on y via £2 7^ 0. 
The insets in Fig. show 7(2:) = 7(2) x 10 4 for the 
values of parameters same as in Fig. ^ (I) V = 0.07, 
r s = 0.2; (II) y = 0.04, r s = 0.6. The function 
7(2) is defined for 2 < 2* and vanishes in the region 
where the plasmon is undamped. In the latter case 
S(z,w) p i = a(z)5(w - Wpi(z)). 

Comparing positions of peaks in constant-2 cross- 
sections of the exact S(z,w) and of the approximate 
S(z, w)so-dam P P i ©, we have observed that the shift of 
the plasmon dispersion due to SO coupling from w p i(z) 
PJl is one or two orders of magnitude smaller than 7(2) 
(unless 7(2) = 0), and therefore can hardly be resolved 
experimentally. 

It is also worthwhile to compare our exact results 
with an approximation commonly used in the litera- 
ture. It is based on the formula H qu » —icr^q 2 / (e 2 u>) 
for small q, where ct w is the ac conductivity (cf., e.g., 
Refs. IB [3). Such approximation takes into ac- 
count only the transitions with zero momentum transfer, 




FIG. 2: The plane (y, r s ) is divided into the domains I and II 
where the plasmon has one or two undamped pieces, respec- 
tively (cf. Fig. 0. Insets show SO-induced plasmon width 
7(2) = 7(2) x 10 4 for the same parameters as in Fig. Q 
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and therefore leads to the kinematic extension of the con- 
ventional particle- hole continuum to a strip y — y 2 <w< 
y + y 2 parallel to the z-axis. The lines w — y ± y 2 co- 
incide with the parabolas w = ±(z ± y) 2 ± (z ± y) at 
z = only, i.e. where the above relation between Tl quJ 
and a u becomes exact. This approximation can be jus- 
tified for very small momenta (e.g., up to z ss 1CP 5 for 
y = 0.01 and up to z » 10~ 7 for y = 0.001). However, 
it fails to describe properly the behavior of the plasmon 
mode at larger z. For example, according to 0] the plas- 
mon spectrum does not cross the border of their damp- 
ing region but rather approaches it exponentially from 



both sides of the strip. This approximation also misses 
the possibility, shown in the left panel of Fig. ^ that 
the plasmon does not emerge undamped for intermedi- 
ate momenta. 

An important test of our results is provided by the 
Kramers-Kronig relations and the sum rules. In partic- 
ular, we have checked numerically that our expressions 
for e(z,w) — 1 and l/e(z,w) — 1 both satisfy the zero- 
frequency Kramers-Kronig relations, which link E2{z,w) 
and S(z 1 w) to the static dielectric function £i(z,0) — 
1 — {r s /z) ■ (2irlli(z, 0)/m*). Our exact result for zero- 
frequency polarization operator is 



27T 7T 

— IIiO,0) = 2 9(1 -y- z) + Q(y- \z - 1|) (1 + - sintp) - 2 <d(z - 1) arccoshz cos if> 



Q(z — (1 + vy)) ( 1 + inp v simp — cos-0„ — 2 cos ip In 



1 + 2 sin(^>„ — vip) \ 
2\/2icos \tj) v cos hip I 



(5) 



where aiatp = y/z < 1 and sini/v = (1 + vy)jz < 1 
are defined in the corresponding ranges of z restricted 
by the unit-step function 0. It is assumed that y < 1/2 
(or 2kfi < hp). For the values z < 1 — y the expres- 
sion iJSJl yields — Hi(z,Q) = m*/n, the density of states 
in 2DEG. In deriving IIi(z, 0) we first perform the mo- 
mentum integration in and only then take the limit 
6 — ► 0. Once this sequence is reversed, there would arise 
the extra term — 7rO(y — z)\J (y/z) 2 — 1 in the rhs of J^J, 
which produces an anomaly at z = y (cf. Ref. |lH)- 

The effect of SO-induced plasmon broadening de- 
scribed above can be observed experimentally. For in- 
stance, one can directly measure S(z, w) by means of in- 
elastic light (Raman) scattering (see, e.g., |l7|). Keeping 
z at a fixed value and varying y and r s , one should ob- 
serve different patterns with either damped or undamped 
plasmon, like those shown in the insets to Fig. ^ where 
the sections of S(z, w) at constant z = 0.1 are presented. 
In order to estimate realistic parameters we focus on 
InAs-based 2DEG where the strength of Rashba cou- 
pling an has quite large values up to ~ 3 x 10 _11 eVm 
|18j and prevails over the Dresselhaus term |l9|. Most 
of experiments deal with the 2DEG's densities rangin 
from n = 0.7 x 10 16 m" 2 13 to n = 2.4 x 10 16 m" 2 " 
Taking Soo ~ 12 and m* w 0.03m e we obtain the range 
of y from 0.04 to 0.075 and the range of r s from 0.38 
to 0.3. These values belong to the domain I in Fig. [21 
where the SO-damping of the plasmon is especially im- 
portant. Choosing the value of 7 ~ 10~ 4 , we establish 
that 2t p i = fi,/(47£i?) is of the order of 10 ps. This can, 
in principle, be resolved experimentally, since the Raman 
measurements done in II- VI quantum wells |2l) give a fi- 
nite lifetime of the plasmon mode with the typical value 



2t p i <~ 0.3 ps. 

Plasmon broadening may be also caused by thermal 
effects. We use the results for 2DEG at finite tempera- 
ture [2^ in order to estimate the relative contributions 
of SO-induced and thermal damping. The former domi- 
nates over the latter below some characteristic tempera- 
ture which is about 90 K for the parameters of Ref. . 
This estimate is made at constant z = 0.1. Note that the 
thermal damping strongly depends on the values of r s 
and temperature whereas the SO-induced damping has 
quite a weak r s -dependence (e.g., compare the insets in 
Fig. 0). This may provide a guide for an experimental 
separation of these two sources of broadening. 

Our result for the dynamic structure factor indicates 
that the other collective mode, longitudinal optical (LO) 
phonon, will also experience a pronounced SO-induced 
damping. Its dispersion and lifetime can be obtained 
from the renormalized propagator which is expressed in 
RPA as El 



D(q,uj) 



2ujlo 



uj 2 -uj 2 n -2u} LO Mm 



(6) 



J LO 

where M 2 = ^AuJLoVq, and A — 1 — Eoc/eq. On the basis 
of this expression we can establish that the LO-phonon 
width (normalized by 2k 2 F /m*) equals 

AS(z,w)w 2 e 2 



1lo{z) 



2we{+ Aw^dei/dw) 



(7) 



w=w ph {z) 



where wq — (m*/2kp) u>lo and w p h(z) is the renormal- 
ized phonon spectrum. Like in the case of plasmons, 
we may neglect the dependence of w p h(z) on SO cou- 
pling, and find the phonon spectrum from the equation 



w 



;[1 + A(l/ei — 1)], where E\ is taken at y = 0. 
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The lifetime effects of the LO phonons can be very im- 
portant for a coupled dynamics of carriers and phonons. 
An interesting possibility arises when uj^q ~ 2orA:.f (or 
w o ~ y)- This, for example, holds in InAs, where the 
value ujlo ~ 28 meV gives Wq ~ 0.07. Changing the 
Rashba coupling strength cir by an applied electric field 
one can manipulate the lifetime of an optical phonon 
which would result in a modification of transport prop- 
erties by virtue of the electron-phonon coupling. 

In Fig. we show the function Jlo( z ) — 7lo(z) x 10 4 . 
In our calculations we used the following parameters: 
w = 0.07, r s = 0.4, eoo = 12 and e = 15 (for InAs). 
Choosing the value 7^0 = 10~ 4 we find that for LO 
phonons the SO-induced lifetime 2t lo is of the order of 
10 ps as well. We note that the typical lifetime for the LO 
phonons measured in AlAs and GaAs by the Raman spec- 
troscopy is of the same order |24|. One can observe that 
for y = 0.07 (\wq ~ y\ < y 2 ) the damping is always finite 
(solid line), while for y — 0.06 {\wq — y\ > y 2 ) it becomes 
nonzero only after some value z « 0.007 (dashed line). In 
the two insets we show the relative location of the phonon 
spectrum (solid line) and of the SO-induced damping re- 
gion (dashed lines). We also depict the boundaries of the 
strip w = y ± y 2 (dotted lines), and one can see that 
for y = 0.07 the phonon spectrum leaves it at the value 
z w 0.017 (left inset), while for y — 0.06 the phonon 
spectrum never gets inside the strip (right inset). 

In conclusion, we have calculated the dielectric func- 
tion of 2DEG with Rashba spin-orbit coupling in RPA. 
We have found that the spin-orbit coupling leads to 
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FIG. 3: SO-induced damping of the LO phonon *f LO (z) — 
7lo(z) x 10 4 for w = 0.07, r s = 0.4, £00 = 12 and e = 
15. Solid and dashed lines correspond to the cases y = 0.07 
(\wo — y\ < y 2 ) and y — 0.06 Qwo — y\ > y 2 ), respectively. 
Insets show the location of the phonon spectrum (solid line) 
relative to the SO-induced damping region (dashed lines) and 
to the strip y — y 2 < w < y + y 2 (dotted lines) for y — 0.07 
(left) and y = 0.06 (right). 



a number of new screening properties: the region of 
particle-hole excitations is extended, and the plasmon 
mode experiences additional broadening. The same 
mechanism leads to the generation of the SO-induced 
lifetime for the longitudinal optical phonons. Thus, spin- 
orbit coupling tends to suppress collective excitations. 
At the same time, it does not affect much the position of 
the collective mode dispersions. Comparison of our theo- 
retical estimates for SO-induced lifetimes with the values 
obtained in the recent Raman scattering measurements 
suggests that the effects described here can be observed 
experimentally. 
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